Synthetic gauge fields and homodyne transmission in 
Jaynes-Cummings lattices 



A. Nunnenkamp^ Jens Koch^, and S. M. Girvin^ 

^ Departments of Physics and Applied Physics, Yale University, New Haven, CT 06520, USA 
^ Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208, USA 

E-mail: andreas . nunnenkampgyale . edu 

Abstract. Many-body physics is traditionally concerned with systems of interacting massive 
particles. Recent studies of effective interactions between photons, induced in the circuit QED 
architecture by coupling the microwave field to superconducting qubits, have paved the way 
for photon-based many-body physics. We derive the magnitude and intrinsic signs of photon 
hopping amplitudes in such circuit QED arrays. For a finite, ring-shaped Jaynes-Cummings 
lattice exposed to a synthetic gauge field we show that degeneracies in the single-excitation 
spectrum emerge, which can give rise to strong correlations for the interacting system with 
multiple excitations. We calculate the homodyne transmission for such a device, explain the 
generalization of vacuum Rabi splittings known for the single-site Jaynes-Cummings model, 
and identify fingerprints of interactions beyond the linear response regime. 



PACS numbers: 42.50.Ct, 42.50.Dv, 71.36.H-c 



Submitted to: New J. Phys. 



Synthetic gauge fields and homodyne transmission in Jaynes-Cummings lattices 



2 



1. Introduction 

The notion of a quantum simulator can be traced back to a remark by Richard Feynman [[T| 
which was later made more precise by David Deutsch [2J and Seth Lloyd [31. Since simulating 
quantum many -body systems on a classical computer is difficult due to the exponential scaling 
of the Hilbert space dimension with the number of particles, it is intriguing to try and use 
a well-controllable quantum system with a tunable Hamiltonian to simulate the physics of 
another quantum system of interest. Over the last decade, this paradigm has motivated 
many experiments and theoretical proposals, predominantly for setups with ultracold atoms 
or trapped ions (see for example Ref . [4J for a recent review) . 

More recently, numerous proposals have emphasized the potential of the circuit QED 
architecture ^ to serve as an interesting solid-state quantum simulator based on polaritons 
ll6l. Indeed, the physics and well-established fabrication techniques for coupled systems 
of on-chip microwave resonators and superconducting qubits, offer a number of promising 
properties. These include an immense freedom in engineering lattices with different 
geometries, dimensionalities, and topologies in a bottom-up approach [|71. Manipulation and 
measurement of quantum states at the single-site level are already standard practice in circuit 
QED experiments. Finally, the controlled coupling to local or global microwave drives also 
promises new insight into the nonequilibrium physics of strongly correlated polariton systems. 

According to theoretical studies [[81 IH [TOl [H] [121, these systems provide access to the 
superfluid-Mott insulator phase transition in the setting of the Jaynes-Cummings model, with 
characteristic similarities and differences compared to the more familiar transition in the 
popular Bose-Hubbard model [fT3l . Theoretical predictions for the physics of finite Jaynes- 
Cummings arrays in the presence of drive and decay exist [[14] [151 [161 [171 [HI [l9l US EH, and 
proposals for introducing synthetic gauge fields [|22l l23l may tum Jaynes-Cummings lattices 
into fruitful ground for studying novel quantum phases of strongly-correlated bosons with 
broken time-reversal symmetry. Experimental efforts to extend the demonstration of single- 
site photon blockade [24] to larger arrays are currently underway. 

In this article, we study the minimal circuit QED system in which breaking of time- 
reversal symmetry is observable: a three-site Jaynes-Cummings ring subject to a synthetic 
gauge field, which induces nontrivial phase factors in the photon hopping terms. Our study is 
organized as follows. In Section 2, we review the many -body Hamiltonian for a general circuit 
QED array, and derive magnitude and intrinsic signs of the photon hopping amplitudes. In 
Section 3, we develop the formalism to calculate linear-response reflection and transmission 
amplitudes. These quantities play a central role in all circuit QED experiments to date, 
which routinely probe quantum states by homodyne measurements. It is to be expected that 
transmission spectra will continue to serve as a fundamental probe in larger Jaynes-Cummings 
lattices, and the first experiments are currently underway^. 

Focussing on the finite Jaynes-Cummings ring, we show in Section 4 the presence 
of degeneracies in the single-particle spectrum for specific values of the gauge-invariant 
phase sum. For commensurate polariton occupation, even weak interactions lead to strong 

I Andrew Houck, private communication (201 1). 
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Figure 1. (a) Jaynes-Cummings ring lattice: photons hop between resonators with amplitude 
t. Each resonator is coupled locally to a superconducting qubit with strength g. (b) Ring 
lattice for ultracold atoms: bosonic atoms can tunnel between adjacent minima of the trapping 
potential with amplitude t and feel an on-site interaction of strength U. 

correlations and generation of highly-entangled Schrodinger cat states. In Section 5, we 
calculate the transmission through such a device. For weak drive strengths, we confirm 
transmission peaks at frequencies corresponding to the eigenmodes within the 1 -excitation 
manifold. Beyond this linear response regime, polariton interactions become important, and 
numerical results show features similar to the nonlinear vacuum Rabi splitting of the single- 
site Jaynes-Cummings model. 

2. Model of circuit QED lattices 

In the simplest setting possible, circuit QED arrays comprise a bosonic lattice with two species 
of lattice sites, cf. Fig. 1(a). The first species, qubit sites, consists of two-level systems 
with transition energy e, realized by suitable superconducting qubits. The second species, 
resonator sites, describes harmonic oscillators with characteristic energy u, implemented as a 
fixed photon mode inside on-chip resonators (e.g. coplanar waveguides). 

Coupling between sites in such a circuit QED lattice allows for photon hopping and 
conversion between photon and qubit excitations, which in part mimics the hopping of bosonic 
particles in the Bose-Hubbard model [fT3l and its multi-species generalizations [l25l . The 
detailed mechanisms for hopping in circuit QED lattices, however, are more diverse: hopping 
t between two resonator sites corresponds to the transfer of individual photons. As opposed 
to quantum tunneling, these processes are possible due to classical EM evanescent fields 
and have a clear correspondence at the classical level. For on-chip microwave resonators, 
this coupling is commonly realized by capacitors interconnecting two or three resonators. 
Hopping between resonator and qubit sites g is of the Jaynes-Cummings type, and enables the 
interconversion of qubit excitations and photons. Direct coupling between qubit sites should 
be feasible for closely spaced qubits, but, as of now, has not been demonstrated in the circuit 
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QED architecture. 

The common situation of small hopping strengths (t, g <^ uj, e) affords the rotating- wave 
approximation, and results in conservation of the total excitation number N = J2j + 
o'j'crj , corresponding to the particle number in Bose-Hubbard type models. In total, the 
Hamiltonian of a circuit QED lattice assumes the form 

j j Ud') i 

We note that disorder in the system parameters (beyond the scope of this paper), may need to 
be included in Eq. ([T]), once it is characterized experimentally. A lattice topology of particular 
interest, is that of a finite-size ring, shown in Fig. [11 a). Its analog for ultracold atoms, see 
Fig. Eb), has been discussed in the literature [|26l |271 |28l |29l|, and (if loaded with bosonic 
atoms) shares some features with the system to be discussed here. 

The option to study broken time-reversal symmetry in circuit QED lattices is based on 
synthetic gauge fields which give rise to complex phases e^^a' in the photon hopping terms. 
Making photons susceptible to a magnetic field is possible, for example, by substituting inter- 
resonator coupling capacitances with passive superconducting circuits [22J. Each such circuit 
consists of a superconducting loop interrupted by three (or more) Josephson junctions, in 
which time-reversal symmetry can be broken by an external magnetic flux. Photon hopping 
takes place via virtual excitations of this coupler circuit, and thus transfers the symmetry 
breaking to the photonic system. Alternatively, active non-reciprocal devices may be used, as 
recently pointed out in an interesting proposal by Kamal et al ll23l . 

2.1. Photon hopping: complex phases and signs 

Individually, the phases e*^Jj' associated with hopping terms are not gauge invariant. Gauge 
transformations of the type dj — )■ djC^^J lead to a fully equivalent Hamiltonian with altered 
phase factors e*^jj'+*^j'~*^^ The relevant quantities are gauge-invariant phase sums, 

^s(c) =g e,,., (2) 

C[j3'] 

around plaquette paths C[jj'] in the resonator sublattice. If no such plaquettes exist (as is the 
case, e.g., in an open Id resonator chain), all phases can be gauged away immediately and the 
hopping terms may be rewritten with positive, real-valued amplitudes, i.e. + 
H.c). For clarity, and to make relevant signs explicit, we will assume t to be real-valued and 
positive throughout; i.e. phases and minus signs will not be absorbed into t. 

As long as time-reversal symmetry is intact, 6y.{C) is an integer multiple of tt on all 
plaquettes, and hopping amplitudes can be chosen strictly real-valued. This, however, leaves 
open an important question about the signs of photon hopping terms. Such signs can arise 
from coupling A/2 resonator modes to each other, or, in general, any half-integer mode with 
L = (n + 1/2)A (n G N). As illustrated in Fig.|2l it is essentially the sign change of the mode 
function for half-wavelength resonances which introduces a gauge-invariant phase 9^ = ir for 
each plaquette with an odd number of resonators. 
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Figure 2. Photon hopping in resonator arrays. Photon lattices, built from quasi- Id microwave 
resonators, admit uniform photon hopping due to capacitive coupling between adjacent ends 
of two or three resonators. When coupling full-wavelength modes (resonator length L matches 
an integer multiple of the wavelength A), photon hopping has a positive sign. For half- 
wavelength modes, L = (n + 1/2) A, the mode function has different signs at opposite ends 
of each resonator, and photon hopping terms can acquire a gauge-invariant phase of tt. It 
is important to note that the Id nature of resonators allows two distinct configurations for 
coupling 3 resonators to each other, see (a) and (b). Gauge invariant tt phases only emerge 
in ring-type coupling (b), and only for rings with an odd number of sites. Rings with even 
numbers of sites, see (c), do not acquire a vr phase. 



More specifically, photon hopping between two adjacent resonators j and f has the 
amplitude 

t e'^ii' = ^CcU (pj (Xend) 0j' (^^end) > (3) 

where Cc is the coupling capacitance and (pjix) denotes the classical mode function of 
resonator j. For full- wavelength modes, the signs of the mode functions can be chosen to 
be positive at each coupling point such that hopping amplitudes are positive. Alternatively, 
for any different sign choice of each mode function (pj, one can always perform a gauge 
transformation that renders all hopping amplitudes real and positive. 

By contrast, half- wavelength modes are associated with an altemating sign, 0j (xend,i) = 
—(pj{xcnd,2), between opposite ends of each resonator, see Fig.|2] The additional signs lead to 
changes in the energy spectrum or band structure if (and only if) the lattice contains at least 
one plaquette with an odd number of sites on its boundary. For such odd-number plaquettes, 
the gauge-invariant phase sum is ^s(mod 27r) = vr. 

Derivation of hopping amplitudes. — We extend the treatment of a single, capacitively 
coupled resonator [|22l to that appropriate for a full circuit QED lattice, and derive photon 
hopping amplitudes and their signs in a circuit QED lattice. Based on the usual circuit 
quantization scheme [|30ll . the Lagrangian for the resonator sublattice is obtained from the 
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LC decomposition of the transmission line resonators and a subsequent continuum limit: 



.(4) 

ends 



Here, $j (x, t) is the voltage profile in resonator j, and i, c denote inductance and capacitance 
per unit length. The sum in the last term runs over nearest-neighbor resonators, Cc is the 
coupling capacitance between each adjacent resonator pair, and the amplitudes are evaluated 
at the corresponding resonator ends connected to Cc- Note that the sign of the coupling term 
stems from the expansion of the capacitive energy |Cc($j — $j')^lends; diagonal terms of the 
former expression appear in the Lagrangian (@J, and also affect the boundary conditions for 
the fields $j [1221 . 

Employing these boundary conditions at the two ports of every resonator, the resonator 
Lagrangian is diagonalized in terms of internal eigenmodes, $j(a;, t) = Cj^(t)e^'^''^(l)j^{x), 
and takes the form 

^ = \ E(C?M - < O - Cc V E <PU^)<PrA^') , i.Ov- (5) 



ends 



In the typical situation of weak coupling between resonators, Cc ^ Cres = cL,'we may focus 
on a single mode, and thus drop mode indices in the following. 

To switch to the Hamiltonian picture, it is convenient to cast the Lagrangian into matrix 

notation, C = ■ {1 — T)C — |w^C^' where Jjj' = 6{jj')Cc(f)j{x)(j)j'{x')\cnds captures the 
coupling between resonators, and denotes the adjacency matrix of the resonator lattice. 
Legendre transform results in the Hamilton function 

n = -n-il-T)-'n--cu''C'' ~ lyf^ - la;^^"'^ + Ivf ■ Tvf, (6) 
2^^ 2 2 2 2 

where the vector vr collects the conjugate momenta hj = dC/dQ- The last approximation in 
Eq. ^ retains the leading order contribution from the resonator coupling. Finally, canonical 
quantization replaces (j — )• {dj + aj) / ^/2u, nj — )■ {ioj — iaj)^fujj2, and in conjunction with 
the rotating-wave approximation, yields the resonator lattice Hamiltonian 

H = ^Y. + ^ H (±)iJ'(«5%' + H-c.)- (7) 
i (i,i'> 

Expressed in proper frequency units, the photon hopping due to direct capacitive coupling is 
therefore given by 

1 



t/27r = -vCc\(l)j{x)(t)f{x')\ 



2 



, (8) 



ends C'rcs 

where u = uj/2n denotes the resonator frequency. As indicated in the discussion above, 
the signs for photon hopping are dictated by the product of mode functions involved in the 
hopping; specifically 

(±)jj/ = sign . (9) 

ends 

This result shows that, up to irrelevant sign changes induced by gauge transformations, photon 
hopping is strictly positive for full-wavelength modes. By contrast, for half- wavelength 
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modes, alternating signs of the mode functions can indeed lead to observable changes. In 
future experiments, this may offer a convenient method for switching between photon systems 
with and without frustrated hopping, by a mere doubling of the microwave drive frequency. 



3. Calculation of steady-state transmission rates 

Experimentally, transmission and reflection measurements are by far the simplest probes for 
circuit QED arrays, and are performed routinely in small circuit QED systems with one or two 
resonators [IBTI l32l l33l . Conceptually, such measurements are implemented by coupling one 
resonator (or several) to an external transmission line, through which microwave photons can 
be injected, and subsequently detected in heterodyne or homodyne. For theory, however, the 
analysis of even this simple setting is challenging for larger numbers of resonators and qubits, 
since it generally involves both strong correlations, and the complexity of non-equilibrium 
physics for a driven, damped quantum system. 

Assuming that the coupling to the external transmission lines (and any additional baths) 
is weak, the system may be described by a Lindblad master equation of the form 

'q = -i [h', ^] + ^ P[a,]^ + 7 5^ I^[<t7]^ + ^ 5Z (10) 

i j j 

Here, the last three terms involve the damping superoperator T>[L]g = LgL^ — {L^L, q] /2 
[[34l and account for photon loss, qubit relaxation and pure dephasing of qubits. For simplicity, 
we will neglect pure dephasing (7^ = 0) and assume decoherence with uniform rates induced 
by coupling to separate baths. Extensions to shared baths and varying decoherence rates are 
possible, but can add further complexity to the analysis. 

The unitary evolution is captured by the Liouville term, and expressed in terms of the 
system Hamiltonian in the frame co-rotating with the drive, i.e. 

H' = UHU^ - tlrU^ = Aat -^dY, -'^'iYl ^t^J + ^d|t=0 = ^0 + H^\t=o- (1 1) 

j j 

Here, the unitary is U = e*'^''* ^ and removes the time-dependence of the drive Hamiltonian 

Hd = J2 (^J- 4^"'"'* + H-c.) (12) 
j 

for a single drive frequency cod, and H = Hi^t + is the Hamiltonian including the drive. 

Under steady-state conditions, ^ = 0, the master equation (fTOl ) forms a linear system 
of equations for the matrix elements of g. Its solution, in principle, enables the prediction 
of transmission and reflection signals by calculating (a) = ti(ag), whose real and imaginary 
part correspond to the two quadratures of the homodyne signal in transmission. It is important 
to stress, however, that the brute-force numerical solution of this equation represents a 
formidable challenge even for arrays with only a handful of resonators and qubits: the number 
of independent real variables to be determined from ^ = is given by (A + i)2A^r22A^q _ 
where Nr and Ng denote the number of resonators and qubits in the lattice, and A is the cutoff 
in photon number. As an example: a lattice consisting of 3 resonators and 3 qubits leads to 



Synthetic gauge fields and homodyne transmission in Jaynes-Cummings lattices 



8 



a linear system of equations with 10^ — 1 variables when admitting up to 4 photons in each 
resonator. A more appropriate cutoff scheme, not based on a truncation of the photon Hilbert 
space, will be described in Section 5. 

In light of this scaling, it is useful to note that the direct calculation of the full reduced 
steady-state density matrix can be avoided, and instead, only those elements of q be evaluated 
that indeed contribute to the expectation value (a). To see this, consider the vanishing time 
derivative of the steady-state expectation value, and re-express it via the master equation: 

d ' / \ Ki 

= ^^^j) = ti{gaj) = -i {^[aj,HQ]j - iVLj - -{aj). (13) 

Not surprisingly, this equation for (aj) in general does not close after explicit evaluation of 
the commutator. For interacting systems it instead produces an infinite hierarchy of equations, 
involving an ever increasing number of expectation values of additional operators, which is 
analogous to situation of the equations of motion method for the system's Green's functions. 
We review the simple and exact solution in linear response first, and discuss numerical results 
beyond this in Section 5. 

Interaction in circuit QED lattices is readily understood as a hard-core repulsion between 
bosonic excitations on qubit sites, which is so strong that levels beyond the ground and first 
excited state are pushed to inaccessibly high energies. Formally, the Jaynes-Cummings lattice 
thus corresponds to a two-species Bose-Hubbard model with Uq = on resonator sites and 
f/i — 7- oo on 2-level sites ll35l . This interaction, of course, is only effective in the presence of 
multiple "particles" (i.e. polariton excitations). Hence, it may be ignored in the linear response 
regime where weak driving does not lead to multiple excitations. In this limit, all qubit sites 
can be replaced by a second set of harmonic oscillator sites, aj — > bj, and the resulting lattice 
Hamiltonian, 

Kt = ^ d]dj + A,J2 + t 5^(e'^^^'aja,v + H.c.) + g ^^(aJS,- + H.c), (14) 

i i UJ') J 

describes a simple system of coupled oscillators. 

Using A„ = v — cudto denote the detuning of the drive from frequency v, the commutator 
required for the calculation of {dj), Eq. (fT3l) . produces 

/iiat]) = A^(a,)+t Yl e''^^'{d,,)+g{b,), (15) 

j'enn(i) 

where nn(j) is the set of sites that are nearest neighbor to site j. Together with the analogous 
relations for d{bj) /dt, one obtains the matrix equation 



9 A, -a 





(16) 



Here, we used the following notation: homodyne amplitudes for the two kinds of sites 
are lumped into the vectors a = ((ai), . . . , (cin)) and (5 = . . . , (&«'))• Likewise, 

VL = (rii, . . . , iln) collects the drive strength for each resonator (dj sites). Finally, T is the 
resonator coupling matrix defined by Tjj/ = te^^n' if j' G nn(j), and = otherwise. 
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From Eq. (fT6] ). one obtains the exact solution for the resonator homodyne signals as 

« = -(A. + T-^;^-zf)"a (17) 

For the physical interpretation of this expression, it is appropriate to switch to the eigenmode 
basis of the resonator subsystem, U(A^ + T)U^ = diag(Ax^), where Ax^, denotes the drive 
detuning from eigenmode frequency T^, and is the matrix of eigenmode amplitudes for the 
resonators; specifically, U*^ = {£u)j gives the relative amplitude of eigenmode u in resonator 
j. Using this basis yields the following expression for the complex homodyne amplitudes 

(%) = - E (At. - k:^:^ - ) ' U., Q,. (18) 

The general structure of this expression can be explained as follows. Coherent driving of 
individual resonators (with amplitudes given by the components of f2) results in an effective 
drive Qr = Uil of the individual resonator eigenmodes. The complex amplitudes induced 
in each eigenmode u depend on the detuning of the drive frequency from photon modes 
and qubits as well as the dissipative rates for photon loss and qubit relaxation. Finally, the 
amplitude measured in resonator j is the weighted sum over all eigenmode amplitudes, with 
weight factors U* accounting for the overlap between resonator j and eigenmode u. 



3.1. Resonator array (g = 0) 

For g = 0, transmission only probes the array of coupled resonators. As expected, a drive on 
a single resonator k results in a transmitted amplitude Tj 1^(10 d) = \ (%) | in resonator j, which 
becomes maximal whenever Ud matches one of the eigenfrequencies T^,. The maximum 
amplitude values (peak heights in the homodyne amplitude) are given by 



At , - iK/2 



(19) 



where fi^^-* = |U*^Ui,fcf2fc| denotes the effective drive strength on resonator j, mediated by 
eigenmode u, and generated by a drive on resonator k. The final approximation of Eq. 
(fT9l ). is valid for non-degenerate and well- separated eigenmodes, i.e. \T^i — t ^ k. 

The generalization to degenerate but otherwise well-separated eigenspaces is achieved by 
including an additional sum over those modes degenerate with u. 



3.2. Jaynes-Cummings lattice (g 0) 

Once the qubit subsystem is coupled to the resonator array, the single-excitation spectrum 
and corresponding transmission resonances are modified. Focussing again on sufficient peak 
separation (g, t 3> 7), the peak frequencies can be approximated by the poles of Eq. (fTSl) 
when setting k = 7 = 0. This yields 

e.,± = i(T. + e) ± i(T, - e)v/l + V/(T,-e)2 (20) 

as the predicted resonance frequencies for the transmitted amplitude. It is simple to verify 
that the frequencies Q^^± are identical to the eigenfrequencies of the full coupled system of 
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the resonators dj and the harmonic modes bj. The corresponding heights of resonances in the 
homodyne transmission amplitude are 

The Jaynes-Cummings physics of resonant vacuum Rabi splitting, and off-resonant dispersive 
shifts evidently carries over to the single-excitation physics of the Jaynes-Cummings lattice. 
The one-excitation spectrum, described by the set of frequencies {Qu,±}, crucially depends on 
the detunings of the qubit frequency e from each resonator eigenmode Tj^. For each individual 
eigenmode u, detuning can be large, |e — T;^| ^ g, and render the interaction dispersive, or 
be near resonant, |e — T,^| <^ g, and lead to strong photon-qubit hybridization. 

Dispersive regime. — The dispersive coupling between qubits and resonator mode u at 
large detuning is reflected in small frequency shifts of the qubit frequency e and the eigenmode 
T^. Expansion of Eq. (|20l) to leading order in g/{e — T^) yields the dispersive shifts familiar 
from the Jaynes-Cummings model, 

Qphoton ^ _ Qqubit^^^^!^^ (22) 

which imply level repulsion between photon and qubit modes. The corresponding height of 
the photon resonance is only slightly reduced from its g = value, whereas the additional 
qubit resonance is strongly suppressed: 

TAQr'') = ^-j^^^^ (23) 

Resonant regime. — Whenever the qubit frequency is tuned into near degeneracy with 
one of the photon eigenmodes, their coupling results in hybridized eigenstates and the typical 
vacuum Rabi splitting. Expansion of equations (|20l) and (|2T1) to leading order in (e — T^)/g 
gives 



1-2 ^ 



Q,^± = T,±g±{e-T,f/g, T,k{Q.,±) ^ 



g{l + k) 



(24) 



7 + K 

for the peak frequencies and heights. 

It is interesting to note that this linearization approach, valid for circuit QED arrays at 
sufficiently small drive powers, may also provide insight beyond the linear response regime 
for multi-level systems with small anharmonicity. The transmon qubit Il32l l36l provides one 
example for such a system, where under appropriate conditions the anharmonicity may be 
treated perturbatively [15J. 



4. Strongly-correlated states in a three-site Jaynes-Cummings ring lattice 

The discussion in Sections 2 and 3 has been general. In the following, we will focus on the 
concrete example of a three-site Jaynes-Cummings ring lattice, see Fig. (2 a). The system 
consists of three resonators coupled via a Josephson ring circuit [[22|. and each resonator 
interacts with its local superconducting qubit. It is the simplest system for which a gauge field 
can lead to measurable effects as the nonzero phase sum around the loop cannot be gauged 
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Figure 3. Three-site Jaynes-Cummings ring lattice, (a) Three resonators are coupled via a 
Josephson ring circuit, and each resonator interacts with its local superconducting qubit. In 
Section 5 we will consider the case where one of the resonators is driven by a microwave 
tone, (b) Spectrum in the 1 -excitation manifold for g = Q (solid) and g/t = 0.2 (dashed). 
For = 7r/3 the highest-energy state is doubly degenerate, (c) Spectrum in the 3-excitation 
manifold. The degeneracy in the 1 -excitation manifold at 6* = 7r/3 induces a 4-fold degeneracy 
in the 3-excitation manifold which is lifted for g/t = 0.5. (d) Overlap of the highest-energy 
3-excitation state with |3, 0,0)| (red solid) and |0, 3,0)| (green dashed) as a 

function of phase twist 9 for g/t = 0.5. Around = 7r/3 the state has large overlap with the 
superposition states \'ip+) in Eq. (|29] |. The insets in (c) and (d) show zoom-outs of the same 
graphs. In all cases, the resonators and qubits are in resonance with each other e = uj. 



away. In this section we show how strongly-correlated states can arise even in this simple 
setup. 

Following Section 2, the Hamiltonian of the three-site Jaynes-Cummings ring lattice is 
given by 

3 

-f^o = ^ l^wo-jaj + ^crjcfj + 9{^]^J + H.c.) 
i=i 

where 9 = 9^/3 denotes the phase twist and 9^: is the gauge-invariant phase sum. We note 
that a similar Hamiltonian arises for bosons in rotating ring lattices Il27ll28ll29ll . 

Since Hq (|25|) conserves the total excitation number N = Yl^j=i + ^t^I^' 

energy spectrum and eigenfunctions can be studied for each separately. In the subspace 



e^^lalai + 0,30,2 + ojos) + H.c. 



(25) 



Synthetic gauge fields and homodyne transmission in Jaynes-Cummings lattices 



12 



with one excitation {N = 1) we can regard the qubits as a second class of harmonic oscillators: 
cTj — 7- bj. Since the system has translational symmetry it is advantageous to introduce quasi- 
momentum operators aj = ^ ^21=0 e^'^'^^^^^^a.k and hj = J2l=o e^'^'^^^'-'^^hk- These reduce 
the Hamiltonian to the form 

i^r^t(4si)(;-)(i:) 

with the eigenmode frequencies = u;+2tcos + 6^ of the coupled resonator sublattice, 
and eigenvalues 

Q± = ± ^ ^. (27) 

Fig. 12b) shows the 1-excitation energy spectrum as a function of the phase twist Q. For 
(7 = 0, the spectrum consists of the resonator modes and the three uncoupled qubits at 
e = cj. For (7 7^ several level crossings between qubit and photon excitations turn into 
avoided crossings, but at the critical phase twist 6* = 0, the ground state (and for 6* = vr/S, the 
highest-energy state) in the 1-excitation manifold remains doubly degenerate, i.e. 6^^ = 6^ 
for 6* = and 6([ = 0^ for 6* = 7r/3. Accordingly, even without synthetic gauge field, = 0, 
the plaquette is frustrated. The degeneracy in the 1-excitation spectrum in turn leads to strong 
correlations among the photons even for weak interaction strength, as we will show next. 

For (7 = each twofold degeneracy in the 1-excitation spectrum induces an (A^ + l)-fold 
degeneracy in the iV-excitation spectrum. Focussing on the critical phase twist Q = n/3, this 
subspace is spanned by \nk=o = n,nk=i = N — n, nfc=_i = 0)| i, i, i) with < n < A^. 
Here, Uk denotes the number of excitations with quasi-momentum number k. 

If (and only if) the number of lattice sites L and the number of excitations A^ are 
commensurate, these degeneracies are lifted by the Jaynes-Cummings term for g ^ 0. To 
understand the induced correlations, we rewrite the Hamiltonian (l25l) as 

2 3 / \ 

Ho = ^ Tfca^afc + e ^ a^aj + — y^^J^^l + H.c. J . (28) 

fc=0 j=l k,j V V'3 / 

The last term lifts the degeneracy between | A^, 0, 0) | |, |, |) and |0, A^, 0) | |, |, \) in a 2A^- 

order process, and the new highest-energy eigenstates are the cat-like superposition states 

|Ar,0,0)±|0,Ar,0) , 
w±) = (29) 

Figures |3lc) and (d) show the spectrum in the 3-excitation manifold {g/t = 0.5, e = to) 
as a function of phase twist, and the corresponding overlaps of the highest-energy eigenstate 
with the states |3,0,0)| and |0,3,0)| 1,4,4). At phase twist < 6* < tt/S the highest- 

energy state has a large overlap with the state |3, 0, 0) | |, |, |). Close to the critical phase twist 
9 = 7r/3, four energy levels approach each other. This is the ( A^ + l)-fold degenerate subspace 
in the g = system. For g ^ the degeneracy is lifted and the two highest-energy states have 
a large overlap with the cat-like superpositions |V^±). For phase twists n/3 <^ 6 < 27r/3, the 
highest-energy state has a large overlap with the state |0, 3, 0) | J,, 4, \). 

We see that synthetic gauge fields can lead to degeneracies in the 1-excitation spectrum 
which may induce strong correlations in the interacting many-excitation system. 
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Figure 4. Transmission in linear response, (a), (b), (d) and (e) show the transmission at port 
2, Ka™*)^' ™ response to driving at port 1, |(a'")p: in (a) for g = as a function of drive 
detuning for A^^ for = (solid) and 9 = tt/Q (dashed) and in (d) for g//t = 1 as well as 
5/t = 1 (solid) and 5/t = 2 (dashed). Transmission for g = as a function of phase twist 
and drive detuning is shown in (b) and for g/ n = 1 as a function of detuning 5 and drive 
detuning A^ in (e). For comparison, we plot the spectrum of the 1 -excitation manifold as a 
function of phase twist 6 for 5 = g = in (c) and as a function of detuning 5 for g/t = 1/3 
and 9 = in (f). The other parameters are t/K — 3 and j/k = 0.1. In (c) and (f) resonator- like 
and qubit-like eigenstates are indicated by solid and dashed lines, respectively. 



5. Transmission of a three-site Jaynes-Cummings ring lattice 



To establish contact with experimental measurements of homodyne transmission, we apply the 
general discussion of Section 3 to the finite Jaynes-Cummings ring. As indicated in Fig. 12 a), 
we consider driving a single resonator with a microwave tone of strength Vl. All resonators 
and qubits are subject to dissipation leading to line widths k and 7, respectively. 

As mentioned above, the eigenmodes of the resonator array are characterized by the 
quasi-momentum number k, and the amplitude of eigenmode k in resonator j is U^.^ = 



0,OVLt 



'2'n-ikj / L 



/vL, independent of phase twist 6. Using the input-output relations 
+ y/KjCij and y/Kjia"^) = i^j [37l, we obtain the reflection and transmission 
rates from the expectation values of the photon annihilation operators, {aj). In reflection, the 
drive interferes with the reflected wave from the cavity, i.e. 



(Si) 



iVt/ K 



and 





2 




2 


1 






Vl/k 







(30) 



In the absence of coupling between resonator and qubits {g = 0) or if qubit dissipation is 
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negligible (7 = 0), we have = (n™) expressing the conservation of the photon flux 

^in/out ^ ^^t) in/out in/out _ 

We first consider the case g = where qubits and resonators are decoupled. Fig. St a) 
shows the transmission |(a2)P as a function of detuning between resonator and drive 
frequency A^^. As expected, peaks in transmission correspond to the eigenfrequencies of 
the resonator array. For 6* = the quasi-momentum states A; = ±1 are degenerate leading to 
only two peaks at —A^/t = — 1 and —A^/t = +2. For 9^0 the spectrum is in general non- 
degenerate and the three peaks mark the three quasi-momentum modes of the ring lattice. In 
Fig. Sfb), the dependence of transmission on drive detuning A^^ and phase twist 9 shows that 
the transmission follows the resonator dispersion ek{9), cf. Fig. life). An interesting detail: 
when the spectrum is almost degenerate transmission is very sensitive to changes in the phase 
twist, see, e.g., the situation of —A^^/t = 1 close to 6* = tt/3. For these parameters the circuit, 
in fact, implements an on-chip circulator (different from the parameter choice of Ref. (TT\). 

For g ^ and small drive, linear response theory remains valid and qubit sites can be 
replaced by harmonic oscillators as discussed in Section 3. Results from this treatment are 
shown in Fig. Hfd) where the transmission 1(02)^ is plotted as a function of detuning A^^ 
(g/t = 1 and fixed 9 = 0). For 5 = e — to = t, the qubits are not in resonance with any of the 
eigenmodes of the resonator array. Relative to the non-interacting case g = 0, the coupling 
to the qubits slightly shifts the frequency of the resonator peaks and three additional peaks 
(two of them degenerate at 6* = 0) appear in the spectrum close to — A^ = 5. By contrast, 
for 5/t = 2 the qubits are in resonance with the k = resonator mode and, analogous to the 
vacuum Rabi splitting in the single-site Jaynes-Cummings model, the resonance peak splits. 
In Fig. Hfe) we show the transmission |(a2)P as a function of the detuning 5 between qubit 
and resonator and drive detuning A^^. As the qubit frequency e is varied, qubits are tuned in 
and out of resonance with each of the resonator modes. For comparison. Fig. Stf) depicts the 
spectrum of the system in the 1 -excitation manifold as a function of qubit detuning 5. 

For stronger drive Vl, population beyond the 1 -excitation subspace sets in, and linear 
response theory breaks down. To find the transmission in that case, we compute the relevant 
expectation values via the steady- state solution of the full quantum master equation (fTOl) . 
As mentioned above, even small resonator arrays and moderate photon number cutoffs lead 
to matrix sizes for the superoperator which quickly become impractical. This problem can 
be mitigated slightly by an adaptive truncation scheme for the density matrix elements. 
First, note that the Hamiltonian Hq conserves the number of excitations A^, while the drive 
Hamiltonian Hd induces transitions between subspaces of different excitation numbers. Each 
element of the standard Fock basis is associated with an excitation number N. Now, 
assigning to each element of the density matrix g^^, = (/i | ^| z/) a generalized excitation number 
N defined as the sum of the excitation numbers of the states and a truncation based 
on N is possible where the cutoff is adjusted with increasing drive strength This cutoff 
scheme is more appropriate and fine-grained compared to a maximum photon number cutoff 
in Fock space. 

In Fig. [3 a) and (b) we present the transmission | (02) P as a function of drive detuning A^^ 
for different drive strengths fi/zt (cutoff: N ~ 4). For weak drive, the response is close to the 
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-A„//; -A^./k -A„/k 

Figure 5. Tmnsmisson beyond linear response, (a) Transmission |(a2"*)/('^i") P ^^'^ (b) 
I (02) P as a function of detuning A^. Linear response approximation (thin black solid in (a)) 
and numerical solution to quantum master equation (fTOt with cutoff N — A for 57/ k = 0.1 
(blue solid), H/k — 0.3 (brown solid with dots), fl/n — 0.5 (green dash-dotted), and fl/K — 1 
(red dashed). Other parameters are t/ k = 5, g/ k = 100, S/k = 5, and j ^ 9 = 0. Blue and 
red vertical lines indicate the one-photon and one half of the two-photon transition energies, 
respectively, (c) For various driving strengths, probability Pq (blue). Pi (green), and P2 (red) 
to have zero, one, and two excitations in the system. Inset in (a) illustrates the fact that the 
2-photon transition energy |0) |2) is detuned from twice the 1 -photon transition |0) |1). 

linear-response result. As the drive increases, the transmission is suppressed compared to the 
linear response value. As better seen in Fig. Eb), the peak additionally becomes asymmetric, 
shifts to smaller detunings and finally at strong drive splits into two. In Fig.[5]we have marked 
the 1 -photon transition energies and half of the 2-photon transition energies, revealing that the 
second peak appears at half the 2-photon transition frequency. This is further corroborated 
in Fig. Oc) where the probabilities for zero, one and two excitations are plotted. At the 
frequency of the 1 -photon transition increasingly strong driving leads to a larger population 
of the 1 -excitation sector Pi but the probability for two excitations P2 remains small. At half 
the 2-photon transition frequency the population of the 1 -excitation space Pi remains small, 
but at stronger drive strength the 2-excitation manifold is occupied. 

This behavior resembles the nonlinear vacuum Rabi splitting as observed in a single- 
site Jaynes-Cummings system [|38ll . Since the spectrum is nonlinear, photon blockade [|24l 
prevents population transfer from the 1 -excitation to the 2-excitation subspace due to the 
energy mismatch of Ei — Eq vs. E2 — Ei, see inset of FiglSfa). Specifically, driving the 
system on the 1-photon transition |0) — )■ |1) transfers population from the ground state |0) to 
the first excited state 1 1), but a second excitation with the same energy cannot enter the system. 
On the other hand, driving the system at half the 2-photon transition frequency |0) — t- |2), not 
very effective at low power, allows for two-photon transitions at sufficiently large power. The 
features of the transmission signal can thus be explained naturally as a consequence of the 
nonlinear spectrum of the system. 
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6. Conclusion 

In summary, we have explored the physics of a Jaynes-Cummings ring, the minimal circuit 
QED system manifesting effects of broken time-reversal symmetry. Analogous to bosons in 
rotating ring lattices, degeneracies at specific values of the gauge field and interactions lead 
to strongly-correlated states. As future experiments are likely to employ transmission as one 
probe for such systems, we have calculated the response to drive and decay, generalizing 
phenomena known from the single-site Jaynes-Cummings model, e.g. the vacuum Rabi and 
nonlinear vacuum Rabi splittings. It is clear that a theoretical description of transmission 
experiments with larger lattices and/or higher drive power will require novel simulation 
techniques beyond the explicit calculation of the steady-state density matrix of the master 
equation. Analytical and numerical methods from the field of quantum transport may turn out 
to be adaptable to this novel setting. Other valuable experimental probes beyond transmission, 
such as two-tone spectroscopy and correlation functions, reveal additional information about 
the system, and will be explored in future work. 
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